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We present a complete set of analytical and invariant expressions for the steady-state density 
matrix of atoms in a resonant radiation field with arbitrary intensity and polarization. The field 
drives the closed dipole transition with arbitrary values of the angular momenta J g and J e of 
the ground and excited state. The steady-state density matrix is expressed in terms of spherical 
harmonics of a complex direction given by the field polarization vector. The generalization to the 
| case of broad-band radiation is given. We indicate various applications of these results. 
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An atomic medium driven by a resonant light field represents a prototype problem in atomic physics and non-linear 
optics. At low density the effects of the atom-atom interaction are small, and the central remaining problem is specified 
by a single atom in a resonant radiation field. As is well-known, the basic processes are the absorption and emission 
I . of photons. The three universal conservation laws (energy, linear momentum and angular momentum) correspond to 
_^ three different aspects in these processes. The main exchange of energy between atom and field corresponds to the 
radiative transition between the atomic energy levels. The momentum exchange gives rise to recoil of the atom, which 
is the basis of the mechanical action of light. The main angular momentum exchange arises from the photon spin. 
Its proper description requires consideration of the light polarization, and the degeneracy of atomic energy levels. 
r"| • Obviously, these processes occur simultaneously, and in a correlated fashion. The recoil effect is usually small, due to 
fS || the smallness of the photon momentum (hk) as compared to typical values of the atom momentum (Mv). In contrast, 
! 1 * . the (spin) angular momentum of photons h is of the same order of magnitude as the internal angular momentum of 
£h ' the atomic states. 

A large fraction of theoretical studies of atoms in radiation fields only considers non-degenerate energy eigenstates, 
without taking into account the magnetic degeneracy of energy levels. In the sense of the group of space rotations, this 
approach corresponds to a scalar model of the atom, which accounts for exchange of energy and momentum, but not of 
angular momentum. This model allows to understand many processes arising from the resonant interaction of atoms 
with radiation [1-5]. However, effects of polarization in combination with the magnetic degeneracy of atomic levels 
cannot be ignored in many cases. The problem remains reasonably tractable in the special cases of linear or circular 
polarization. In these cases, there is an obvious choice of the quantization axis, so that the submatrices of the density 
matrix for the ground and the excited level remain diagonal at all times. For arbitrary elliptical polarization the 
situation is appreciably more complex. There are various experimental situations where elliptical polarization is quite 
essential. An important example is the phenomenon of coherent population trapping (CPT) of atoms driven by an 
elliptically polarized light field [6]. Another example is cooling and trapping of atoms in light fields with polarization 
gradients [7] , where the strong correlation between the processes of linear and angular momentum exchange can lead 
to atom temperatures down to the single-photon recoil energy (~ lCT 6 !i). Here continuous spatial variations of the 
polarization are crucial, which, except for special cases, give rise to elliptical polarization. Models with non-degenerate 
atomic states can only describe the Doppler limit of laser cooling (~ 10~ 3 K) [5]. 

A central part of these processes is the resonant interaction of an elliptically polarized light field with a closed 
atomic transition between degenerate energy levels. In this case the light-induced anisotropy of the atomic state 
is long-lived. This enables one to accumulate information on very weak couplings, which allows for high-resolution 
spectroscopy. As noted above, for many cases one can consider the recoil effect as a small perturbation of the order 
of hk/Mv <C 1. In zeroth order the atom has a constant linear momentum. In this case only the exchange of energy 
and angular momentum between atom and field is taken into account, and the state of the atom is fully described 
by the density matrix for the internal states. The remaining field effects (light shifts, field broadening, change of 
population and coherence etc.) are caused by the stimulated and spontaneous transitions, which are described by the 
generalized optical Bloch equations (GOBE) for the internal atomic density matrix. Depending on the light intensity 
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and interaction time, three limiting cases can be distinguished. The first case occurs for short light pulses, when 
the interaction time is so short that relaxation processes can be neglected. Then the interaction of atoms with the 
field has a coherent character, and it can be described by the time-dependent Schrodingcr equation for the atomic 
wavefunction. Typical effects are coherent transient processes, such as Rabi oscillations or photon echoes, which have 
been thoroughly analyzed [1,8], with or without inclusion of the magnetic degeneracy of atomic states [9]. The second 
case occurs when the interaction time is long compared with the spontaneous lifetime, but still short compared with 
the absorption time, which determines the rate of optical pumping. In this case one can use perturbation theory. This 
situation arises for laser beams of low intensity and small diameter, when time-of-flight effects become important. 
Solutions of the optical Bloch equations corresponding to this case have been obtained with the use of irreducible- 
tensor techniques for arbitrary Zeeman and hypcrfine level structures [10]. Finally, when the interaction time is so 
long that perturbation theory becomes inapplicable, one has to find the steady-state solution of the GOBE. This 
situation appears either for slow atoms interacting with light, such as in optical molasses, or in the case of light beams 
with high intensity or large diameter. A common restriction is that only a closed transition between two degenerate 
atomic levels is considered, so that the total population is conserved. 

The GOBE are an essential ingredient of the description of sub-Doppler laser cooling by fields with polarization 
gradients [11]. We here recall just a few representative cases. In the semiclassical theory of laser cooling authors 
concentrated their efforts on the velocity-dependent steady-state density matrix. The presented analytical results 
were, however, restricted to transitions with specific small values of the angular momenta of the states [12]. Berman 
ct al. [13] have formulated the GOBE for arbitrary field polarization and arbitrary atomic energy level structure, 
using the irreducible tensor representation. They demonstrate that the sub-Doppler light forces and the sub-natural 
resonances in non- linear spectroscopy are closely related. Similar results have been obtained in [14], where a general 
relationship between the light force and the non-linear polarizability tensor has been derived. 

For linear or circular polarization, the steady-state solutions of the GOBE have been discussed in various papers 
[15-17] for arbitrary values of J g and J e . For arbitrary polarization, the symmetry is reduced, and the steady state 
for arbitrary J g and J e represents a complex mathematical problem. The number of equations, which is equal to 
the number of elements of the density matrix, amounts to 4(J g + J e + l) 2 . The steady-state density matrix was 
found in analytical form only for transitions involving specific small values of the angular momentum (J g = 0, 1/2, 1) 
[18-22]. Recently we have discussed the steady state of atoms in light fields with arbitrary polarization, for transitions 
with J e — J g = [23,24] or 1 [25]. In these cases, the structure of the solutions looked remarkably different. An 
invariant approach to the general problem, based on the expansion of the density matrix in bipolar harmonics of 
complex directions was developed in [26]. Nasyrov [27] has suggested an alternative approach to the steady-state 
density matrix, using the semiclassical Wigner representation of angular-momentum orientation . This method seems 
especially useful at large angular momentum J ^> 1. The results are in good qualitative agreement with our exact 
solution. 

In the present paper we give unified exact analytical expressions for the steady state of atoms driven by light with 
an arbitrary polarization, for all possible dipole transitions. Radiative relaxation is included in the description, and 
the results are presented in an invariant form. The analysis is based upon the group-theoretical properties of the 
transition dipole. The polarization direction in real space is reflected in the density matrix for the two subspaces 
corresponding to the ground and excited state. The symmetry group for the problem is the group of rotations SU(2), 
which is an important leading principle for the search of the solution, as well as for its presentation in invariant form. 
The general discussion allows us to clarify the similarities in the different cases, which were not obvious at all in the 
separate treatments. 

The nature of the steady state strongly depends on the value of J e — J g . For dipole-allowed transitions, this 
difference can be —1, or 1. For the case that J e = J g = J we should moreover distinguish the cases that J is integer 
or half-integer. This leads to four classes of transitions. 

a) Transitions J g = J — > J e = J — 1. In this case atoms are optically pumped into dark states, where they do not 
couple to the light field. This is the phenomenon of coherent population trapping. These dark states are linear 
superpositions of the Zeeman substates of the ground level, defined as eigenstates of the resonant interaction 
Hamiltonian with zero eigenvalue. Hence, there are no light shifts. For the present class of transitions, there are 
two independent dark states, which span a two-dimensional space. This space depends only on the polarization, 
and it is independent of the field intensity and the detuning. Both the atomic dynamics in the field and the 
steady state depend on the initial state. 

b) Transitions J g = J — > J e = J with J integer. In this case a single dark state exists, and CPT takes place, so 
that in the steady state the atom is in this unique pure state. Therefore, the steady state does not depend on 
the initial conditions, the intensity or the detuning. 

c) Transitions J g = J — > J e = J with J half-integer. For this class no dark state exists, and CPT does not occur. 
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The only exception is the case of circular polarization, where a single dark state does occur. The steady state is 
uniquely defined, but now it depends both on the polarization, the intensity and the detuning. Moreover, it is 
not a pure state. In fact, it has the remarkable property that the excited-state submatrix of the density matrix 
is fully isotropic, which makes the analytical expression for the entire density matrix particularly simple. 

d) Transitions J g = J — » J e = J + 1. For this class of transitions the steady-state solution is unique. There is no 
dark state. The excited-state submatrix is always anisotropic. 

Only in the cases c) and d) does a steady-state excitation exist. In both cases, the anisotropy of the excited state 
and of the optical coherences depends only on the polarization of the driving field. The intensity and the detuning 
enter only as an overall multiplicative factor. Moreover, in both cases we find that the submatrices both for the 
excited and the ground state are even functions of the detuning. In the cases a) and b) only the ground state is 
populated in the steady state, and the excited-state submatrix and the optical coherences vanish. The occurrence of 
dark states and velocity-selective CPT allows to reach cooling below the recoil limit [28], and it has been studied by 
many authors [29] . In the case of J g — 1 an invariant form of the dark states in an elegant vector notation has been 
used for the analysis of CPT in 2D and 3D [30]. Here we extend an invariant approach to all the dark transitions. 

The remainder of the paper is organized as follows. First we discuss the general structure of the generalized optical 
Bloch equations (Sec. II), which leads in Sec. Ill to the definition of a natural basis of states that depend only on 
the light polarization. We separately discuss the cases with (Sec. IV) or without (Sec. V) dark states. In the latter 
case, we discuss the effect of optical pumping on the degree of excitation and the AC Stark shift in Sec. VI. The 
generalization to the situation of broad-band radiation is given in Sec. VII. 

II. FORMULATION OF THE PROBLEM 

We consider a closed atomic transition J g — > J e with a transition frequency luq of an atom at a given position. In 
the present paper we will not consider the translational motion of the atom. This corresponds either to the case of 
very slow atoms or to the case of a travelling plane wave, where the atomic motion at a given velocity leads only to a 
Doppler frequency shift. The transition is driven by a monochromatic radiation field with frequency u> and arbitrary 
polarization e. The time-dependent electric field vector at the position of the atom is given by 

E = E e cxp(-iujt) + c.c. (1) 

where the polarization vector is 

e = E e ^ = E (- 1 ) 9e -^- ( 2 ) 

q=0,±l 9=0, ±1 

Here E is the complex field amplitude, e q = e • e q is the covariant spherical component of the polarization vector 
e, and the spherical basis vectors of polarization are defined by {eo = e z : e±i = =F( e a; ± ie y )/\/2}. Notice that 
e* = (— l) 9 e_(j. The vector e is normalized, so that e* • e = 1 and without loss of generality we assume that its real 
and imaginary parts are orthogonal, which implies that Im(e • e) = 0. Then the two vectors Re e and Im e are the 
axes of the polarization ellipse. 

It is always possible to use a coordinate frame where only two of the components e q are nonzero. There are two 
possibilities. The conventional choice is that the Oz axis is chosen normal to the polarization plane. In this coordinate 
system the vector e is the sum of the two opposite circular unit vectors e±i. If the Ox axis is directed along the 
major semiaxis of the polarization ellipse (see Fig. la), e is written as 

e = e-r cose + iej, sine = — e+i sin (e + tt/4) + e_i cos (e + 7r/4) , 

where the ellipticity angle e can take the values — 7r/4 < e < 7r/4. Obviously, | tan e| is equal to the ratio of the minor 
semiaxis to the major semiaxis and the sign of e determines the helicity. 

Another choice is called the natural coordinate frame, which was introduced in [31]. When we represent the 
polarization ellipse as the intersection of a cylinder with a plane, the natural frame Ox'y'z' is specified by the 
requirement that the axis Oz' is the axis of the cylinder, while the axis Oy' coincides with the axis Oy. The minor 
semiaxis of the ellipse coincides with the radius of the cylinder (Fig. 2a) . Then the polarization e is the superposition 
of a linear component along Oz', and one circular component. The two frames are connected by a rotation along the 
axis Oy over an angle 9 obeying the relation 

cos — | tan e\ . 
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In the natural frame, the polarization vector is specified as 



e = eQ-\/cos(2e) — e' ±1 V2am(e) 



(3) 



where the helicity of the spherical unit vector e' ±1 corresponds to the sign of e. In general there are two possible 
choices for the cylinder, corresponding to opposite signs of the rotation angle 9. Notice that when the polarization e 
is represented in terms of the Stokes vector as a point on the Poincarc sphere, the angle 9 is equal to the polar angle 
of this point [32] . 

The quantum kinetic equation for the density matrix p of the internal state of the atom in the external field (1) 
has the form: 



d 



H ,p 



-d-E(t),p -r{p}. 



(4) 



Here H is the Hamiltonian describing the energy of the two resonant levels of the free atom and d is the dipole 
operator connecting the two levels. The radiative relaxation is described by the operator r{p}. All operators are 
represented as matrices on the Zeeman basis of the ground and excited levels, with states {\J g ,p g )} and excited 
{\J e ,p e )}. The density matrix p can be separated in four matrix blocks, where the matrices p gg and p ee are the 
submatrices for the ground and excited state, and the off-diagonal blocks p eg and p ge describe the optical coherences. 
In the rotating-wave approximation the time dependence of the kinetic equation can be removed by introducing the 
transformed optical coherences as 



p eg = cxp(-iujt)p eg p ge = exp(iujt)p : 



gc ■ 



(5) 



The resulting system of generalized optical Bloch equations (GOBE) can be expressed in the dimensionless dipole 
operator D, which couples the ground state to the excited state. It is specified by the definition of its spherical 
components as 



(6) 



(M) 



so that their matrix elements are equal to the Clebsch-Gordan coefficients C J T el f: e , „. Stimulated transitions are 

Jg fig Lq 

described by the operator V, which is the component of the vector operator D in the polarization direction: 



V = D • e . 



(7) 



Then the GOBE take the form 



d_ 

dt 



-iS 



eg 



V Pgg ~ PeeV 



(8) 



d 7 

ai + 2 +tS 



= in* [yt;5e 



Pgg 



(9) 



dt 



+ 7 p. 



(10) 



-Pgg-l 2^ 

9=0, ±1 



dt 



Dg=i 



ge 



(11) 



with the normalization 



Tr{p s J + Tr{p ee } = l. 



(12) 



Here S = lu — oj eg is the detuning, uj eg = (E e ~ E g )/h is the transition frequency, 7 is the radiation relaxation rate and 
O = Eo(J e \\d\\J g ) /h is the generalized Rabi frequency, expressed in the complex field amplitude Eq and the reduced 
dipole matrix clement (J e ||e?||J a ) that determines the strength of the transition. The summation in cq. (11), which 



4 



describes the feeding of the ground state by spontaneous decay, runs over the three possible independent polarizations 
{eo,e±i} of spontaneous emission. Conservation of the total population of the closed transition is ensured by the 
relation 



9=0, ±1 

When acting on an isotropic excited state, this feeding term is proportional to 

V &d - 2Je + 1 n 

g=0,±l M 



(13) 



(14) 



We introduced IT g and IT e as the projectors on the ground state and the excited state. The dynamical equations (8)- 
(11) represent the generalized optical Bloch equations, which describe transient processes as generalized damped Rabi 
oscillations, optical nutation, free induction decay, etc., as well as optical pumping effects, that lead to an anisotropic 
distribution of atoms over the magnetic sublevels. 

Equations for the steady state are obtained by setting all time derivatives to zero in eqs. (8)-(ll), which gives a 
set of linear equations for the density-matrix elements. By using eqs. (8) and (9), the steady-state optical coherences 
can be directly expressed in the population submatrices p ee and p gg as 



-in 



y/2-iS 



Vpgg - p ee V 



-in* 



Pee Pgg 



rge 1 /2 + i8 

After substitution in (10) and (11), this leads to closed equations for the population submatrices in the form 

jp ee = - 1 S/2{VV\ p ee } + jSVp gg V^ +i5S[VV\p ee ] , 



(15) 



(16) 



-7 £ Dl% e D q = - 7 S/2{V^V,p gg }+ 1 SV% e V-i5S[V^V,p gg \; (17) 

9=0, ±1 



where {,},[,] indicate an anticommutator and a commutator, respectively, and where 

\nf 

7 2 /4 + <5 2 



(18) 



is the saturation parameter, which is proportional to the light intensity and to the global oscillator strength of the 
transition. The left-hand sides of eqs. (16) and (17) describe spontaneous processes, i.e. the radiative damping of 
the excited level and the spontaneous transfer of population and Zeeman coherence from the excited to the ground 
level. The terms in the right-hand sides that are proportional to the optical pumping rate "/S represent light-induced 
loss (with the minus sign) and gain (with the plus sign) of the levels. The commutator terms in the right-hand sides, 
proportional to 5S, describe the AC Stark effect. They contain the light-shift operators in the ground and excited 
level 

£ g = SSV^V , £ e = -5SVV ] , (19) 

which play the role of an effective Hamiltonian for the two levels. 

The steady-state solution of the GOBE corresponds to the limit t — > 00, with t the interaction time. In practice this 
means that t is larger than the largest relaxation time in the internal degrees of freedom. In the case of a degenerate 
ground state at low saturation this largest time is of the order of {^S)~ l , which is the inverse of the rate of optical 
orientation in the ground state. For large saturations S > 1 the largest relaxation time is the excited-state lifetime 
7 _1 . Thus, the conditions for the steady-state regime can be written as 

min{7i, 7^} > 1 . (20) 
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III. NATURAL BASIS OF STATES 



A. Commutation of density matrix and light-shift operators 

In the special case of linear or circular polarization, the Zceman substates \J g ,p, g ) and |J e , p, e ) constitute an obvious 
natural basis of substates in which to express the density matrix. For linear polarization, one chooses the quantization 
axis parallel to the polarization direction, so that the polarization vector e is equal to the spherical unit vector e q 
with q — 0. For circular polarization, the polarization vector e is equal to the spherical unit vector e q with q = ±1, 
provided that the quantization axis is chosen normal to the polarization plane. In both cases, the operator V couples 
each substate \J g ,n g ) to a single excited state |</ e ,/i e ), which /i e = p g (linear polarization) or ji e = p, g ± 1 (circular 
polarization). In this case it can be easily checked from eqs. (8)-(ll) that the equations for the populations do not 
mix with those for the Zceman coherences. Also, eqs. (16) and (17) show that the steady-state solutions p ee and p gg 
are diagonal on the basis of the Zeeman substates | J, p) . Since also the light-shift operators £ g and £ e are diagonal 
on the Zeeman substates, this implies that the steady-state density matrices p gg and p ee commute with the light-shift 
operators £ g and £ e , so that for linear or circular polarization we find 

[p gg ,£ g ] = 0, [p ee ,£ e ]=0. (21) 

Moreover, when all Zeeman coherences are zero initially at t = to, they remain zero for all later times t > t - 

One might be tempted to believe that also for arbitrary elliptical polarization a basis of states can be chosen 
for which populations and coherences do not mix. However, this is not true. It has been shown in ref. [33] that 
spontaneous decay can create coherence between eigenstates of £ g , even if they do not exist initially. In general the 
time-dependent solutions p gg and p ee for arbitrary polarization will not commute with the light-shift operators at all 
times. 

Nevertheless, in this paper we shall prove that for the steady-state solutions for all classes of transitions the 
commutation rules (21) are valid for arbitrary elliptical polarization and for all classes of transitions. The proof is 
rather different for the various classes, so that it is most convenient to give the proof while discussing the expression 
for the steady state for each class separately. An immediate consequence of the commutation rules (21) is that the 
steady-state density matrix is diagonal in the eigenstates of the light-shift operators. This implies also that the last 
terms in eqs. (16) and (17) vanish, so that the steady-state population submatrices depend on the detuning 5 and 
the spontaneous-decay rate 7 only through the saturation parameter 5*. 



B. Eigenbasis of light-shift operators 

We are interested in the eigenstates of the operators V^V and VW . The corresponding eigenvalues are real and 
non-negative, so that we can write 

V^V\{9)i) = >H\{9)i); VV^\{e)j) = \%e)j) (22) 

with Xi real, and with eigenstates \(g)i) in the ground and |(e)j) in the excited level. For given values of J e and J gi 
the eigenstates and eigenvalues are fully determined by the polarization of the driving field, and they do not depend 
on the detuning or the intensity. The states |(e)j) and \(g)i) form the natural bases for the excited and the ground 
level. The operators £ g and £ e are diagonal with diagonal elements 5SXf and — SSXj. Operating with the coupling 
matrices V and V* on the first equation (22) shows that V\{g)i) is eigenstate of VV^ with eigenvalue Xf. Hence we 
may assume that V\{g)i) is proportional to \{e)i). A proper choice of the phases of the eigenstates \{e)i) then leads 
to the expression 

V = ^X i \(e)i)((g)i\. (23) 

i 

Hence each non-zero value of Xi corresponds to a pair of states | (g)i) and \(e)i) that are coupled by V and . In 
addition, the operator VW or VV^ may have eigenvalues zero. The corresponding eigenstates are unaffected by the 
radiation field, and they do not contribute to the coupling operator (23). 

If the commutation rules (21) are true, the steady-state density matrices p gg and p ee are diagonal on these bases. 

The diagonal elements 7rj and ir^ 1 are the stationary populations. Taking the diagonal elements of the equations 
(16) and (17), we obtain the relations 
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j4 e) =-lS\rt e) +lS\M 9) (24) 
j 

for the steady-state populations, with 

Wi,= E \((e)j\D q \(g)i)\ 2 

q=0,±l 

the probabilities of the spontaneous transitions j —* i. These transition probabilities are normalized as ^ i Wy = 1 
for all j. As is seen from (24), if £ e has an eigenvalue equal to zero, then the corresponding eigenstate |(e)j) is not 
populated. Conversely, if one or more eigenvalues Xf of £ g are equal to zero, then a steady state exists where only the 
corresponding eigenstates \ (g)i) of the ground level are populated. For Xj ^ it follows from (24) that the populations 
of the ground- and excited-level substates are related by the equation 

^ (9) = (l + ^M e) - (26) 

From (25), it then follows that one can deduce a closed system of equations for the excited-state populations in the 
form 

£ Wtf7r <«> =7r <«>. (27) 

j 

These relations (27) uniquely determine the populations 7rj e \ apart from normalization. This implies that the steady- 
state density matrix of the excited level p ee depends on the intensity and the detuning only through a normalization 
constant, which is a function of the saturation parameter S. Moreover, they show that the distribution over the 
excited-level substates can be considered as a stationary point of the radiative relaxation operator, in the sense that 
such a distribution is invariant under spontaneous decay to the ground level [34] . 



C. Condition for diagonal steady state 

The conjecture of the commutation relations (21) can be formulated in an invariant form. It is sufficient to assume 
the existence of two Hermitian operators E and G, with E acting on the excited states, and G on the ground states, 
and obeying the identities 

EV = VG , D\ED q = G . (28) 

<?=o,±i 

From the first identity (28) it follows that the operators E and G have an identical diagonal matrix form on the 
natural bases. From this identity and its Hermitian conjugate V^E — GV^ one obtains the commutaton rules 

[E,VV*] = , [G, = . (29) 

Starting from the relations (28), while using the equations (16) and (17), one easily derives that the steady-state 
submatrices p ee and p gg are determined by the relations 

p ee = PSE , V^Vp gg = p (l + SV^V) G , (30) 

with P a normalization constant. These relations are just the operator expression of eq. (26). From (29) and (30) 
it follows immediately that the commutation rules (21) hold, and that the density matrix is diagonal on the natural 
basis. Therefore, the problem of finding expressions for the steady-state density matrix is now reduced to finding 
operators E and G that obey the relations (28). These operators can be assumed to depend only on the polarization 
vector e, and not on the intensity or the detuning. When operators E and G obeying (28) are found, the steady-state 
density matrix is determined by (30), and it is indeed diagonal in the natural basis. 
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IV. DARK STATES 



As recalled in the Introduction, dipolc transitions can be classified into two classes depending on the occurrence 
of coherent population trapping (CPT). For the first group, where CPT occurs, one or more of the eigenvalues of 
the ground-state light-shift operator £ g vanish, so that this operator cannot be inverted. During the optical-pumping 
process atoms are accumulated in the corresponding eigenstates, which are termed dark states, since they do not 
interact with the light field. Then the trivial solution E — G — of the system (28) still determines a normalizable 
steady-state solution of the relations (30) obeying p ee = and V^Vpgg = 0. The ground-state density matrix p gg is 
composed of the ground-level dark states \^^ NC ^), which obey the equation 

^|*( JVCf )) = 0. (31) 

In order to specify the dark states in an invariant manner, we view state vectors as tensors. A state vector |W) in 

the ground level is considered as a tensor Hfj g of the rank J g , with (covariant) components _ Ms specified by the 
expansion 

\*) = J2(-i)-w Jg -, 9 \j g ,» g ) . 

Ms 

Using the Wigner-Eckart theorem, we express the matrix elements of the left-hand side of eq. (31) as 
(J e ,Me|(d • e)\¥ NC ^) = (J e \\d\\Jg) E ° J J'X 



,(JVC) 
Jg — Ms 



<?,Ms 



= <J e ||d||J fl >(-l)-'"{e®*^^} Je _ /le) 

where {. . . <£> . . .} denotes the standard definition of an irreducible tensor product [35], (J e ||c£||J g ) is the reduced 
matrix element of the dipolc moment operator, and as before, e q = e • e q is the covariant spherical component of the 
polarization vector e. The invariant expression of eq. (31) in terms of a tensor product of rank J e reads 

{e®*£ C) } Je =0, ( 32 ) 
with the normalization condition (^f C ^ ■ ^j^ C ^) = 1 . 



A. Transitions J g — J — > J e = J with integer J 

For transitions with integer values of J g = J e = J there is a single dark state for any polarization [6]. In order 
to write an explicit and invariant form of tyy c ' in this case, we introduce the L-fold tensor product of the vector e 
[36,37] 

{e} L = {. . . {{e ® e} 2 ® e} 3 ... ® e} L . (33) 

which are proportional to the spherical harmonics of a complex direction 71lm(g) (Appendix A). Notice that the 
three components {e}i 9 in the case L = 1 coincide with the spherical components of the (possibly complex) vector e. 

The Clebsch-Gordan expansion (A. 3) for a product of two spherical harmonics with the same argument leads to 
the result 

{{e} L {e h } K = C« \ ^ ^ (2L _ i y (2J l 1) ,, {-}k • (34) 

It follows from the symmetry of the Clebsch-Gordan coefficients that C^§ J0 = 0ifL + J — K is odd. Then choosing 
L = 1 and K = J we obtain {e (g) {e}j}j — 0, so that the single dark state as defined by (32) can be specified in 
tensor form as 

^=M{e}j. (35) 
The normalization constant follows from the equality 
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«•>' ■ <•">'> - Ira' 6 ■ e) '" v ■ e ' y ' 2 Pj ($3?) • <*> 

as an example of the sum rule for spherical harmonics (A. 4). Here Pj(x) is the standard notation for Legendre 
polynomials. This leads to the expression 



.(2J-1)!! v ' V( e ' e ) 
iiati* 

harmonics. The steady-state density matrix 



'1/2 

(37) 



In general, the algebraic and transformational properties of the dark state ^>^ c ^ are the same as for spherical 



p 9S H* (JVC) ><* (JVC) l (38) 

obviously commutes with the light-shift operator £ g . 

In the special case of J g = J e = 1, which is the prototype case of CPT [28,30], the dark state is specified by 

^(nc ))= £ (-1)^1,^ = 9). (39) 

9=0, ±1 

It is well-known that the states \ l,q) with angular momentum 1 have the same transformation properties as the three 
spherical unit vectors e q , so that any state vector can be represented by a Cartesian vector. When the state coupled to 
an excited state with angular momentum 1 by the operator V, the vector representing the excited state is represented 
by the vector that is the cross product of the ground-state vector and the polarization vector, since the cross product 
is the only way in which a vector can be formed from two vectors. Now a comparison of (39) with (2) shows that the 
expansion coefficients are identical, so that the state (39) has the polarization vector e as its vector representation. 
This immediately explains why eq. (31) holds for this state (39), since the cross product of a vector with itself vanishes 
[30]. The explicit invariant form (35) of the dark states for integer values of J g = J = J e generalizes the well-known 
result for J = 1. 



Je = J - 1 

For transitions with J e = J g — 1, the CPT-condition (32) takes the form 

{e®* { J NC) }j- 1 =0. (40) 

In this case there is a two-dimensional dark subspace, spanned by two independent dark states [6]. The scheme of 
light-induced transitions in the natural coordinate frame is shown in Fig. 2b. It is explicitly seen that the dark state 
coincides with the outermost Zeeman substate \J,fj, = J). The other linearly independent dark state is determined 
analogously in the frame, connected with the second cylinder. 

First we consider the case that J is integer. An invariant tensorial expression for the dark state is directly obtained 
when we notice that the outermost Zeeman substate is given by the tensor {C} j with C the circular component of 
the polarization vector in the corresponding natural coordinate frame. This vector C is completely specified by the 
requirements that it is normal to the polarization vector. Hence, the solution of (40) is represented by the tensor 

^j NC) = {C}, , (41) 

where 

(C*-C) = l ; (C-C) = ; (e • C) = . (42) 
The two independent (but not orthogonal) solutions are 

C (i,2) = [ex [exe^^fexe-yfe-e) 
V(l-|e-e|2)(l + |e-e|) ' 

and the two corresponding dark states are called vE'j'' and . These states are normalized and linearly independent 
but not orthogonal. We can combine them into two orthogonal states defined by 
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*(±) = J J (44) 

/ 2 (l ±( ^.^))) 



In order to calculate the dot product ■ ^[p^j wc take into account the relationship 

((a), ■ (b W - ^ (v^Hb-TT)'' P., ( V( .^. b) ) - • •>)' ■ 

which holds if either a or b is a circular vector. Hence we obtain 

(,<»■.,») _(o«-.o») J -(l^)' («, 

Next we turn to the case of half-integer values of J. Then we can use the correspondence between circular vectors 
and spinors. The tensor product of a spinor \ (defined as a tensor of rank 1/2) with itself into a tensor of rank 1 is 
always a circular vector, so that 

{x®xh=C, (46) 

with C • C = 0. The plane of this circular vector is normal to the direction of the orientation of the spin vector 
represented by the spinor. Conversely, any circular vector can be represented in the form (46) for some spinor \- Now 
for a given polarization vector e, the two circular vectors (43) correspond to two spinors x' 1 ' 2 ) so that 

} = C [m) , m = 1,2 . (47) 



x (m) ® x (m) 



Since the two dark states <J/^™- ) ( m = 1, 2) are the outermost Zeeman states in the two natural coordinate frames, 
they can be expressed in the form 

* ( 7 l) = {X [m) }j , (48) 
where the tensor {x}j is constructed from 2J spinors \ 

{x}j = {••■{{x®x}i®x} 3/2 -'®x} J • 

The orthonormalization is specified by equations (44) and (45) for both integer and half-integer momenta. 

In the steady state, the excited submatrix p ee disappears, and the ground-state density matrix p gg can be an 
arbitrary density matrix within the two-dimensional subspace spanned by the two dark states 



\*i NC) ) = T,(-^^ 9 \j a ,» 9 ). 

fJ-g 



Obviously, any density matrix within this subspace commutes with the light-shift operator £ g . For any value of J, 
this dark subspace depends only on the polarization vector e. However, the specific steady-state density matrix in 
which an atom will end up can depend upon the initial state as well as on the intensity and the detuning of the light 
field. This case of a transition with J e = J g — 1 is the only case of a dipole-allowed transition in which the steady 
state is not unique. 

V. NO DARK STATES 

A. General form of steady state 

For a transition without dark states, the steady-state solution is unique, as has been proved in Ref. [41]. Then the 
excited level is populated in the steady state. This is the case when the ground-state light-shift operator £ g has no 
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eigenvalues zero, so that the operator V s V acting within the 2 J g + 1 states of the ground level can be inverted. When 
operators E and G exist that obey the relations (28), the steady-state is obtained from (30) in the form 

p ee = /3SE , p gg = (3 [(V^V)- 1 + S\G (49) 

From the commutation rules (29) it follows that the density matrix is diagonal on the natural basis. The optical 
coherences are directly evaluated from eq. (15), with the result 

^=(M t= (^W^- lg - < M > 
The constant (3 follows from the normalization condition (12), and we obtain 

0= * , (51) 
cto + Zbai 

with the invariant expressions for the coefficients 

a t = Tr{E} = Tr{G} (52) 

and 

ao-Trjc^)-^} . (53) 

These coefficients «o and ot\ depend on the polarization only, not on the intensity or the detuning. The steady-state 
density matrix depends in the intensity and the detuning only through the value of the saturation parameter S, 
defined in (18). 

It is noteworthy that the submatrix p ee of the excited level is always proportional to the single operator E. 
This implies that the steady-state anisotropy of the excited level, such as its orientation and its alignment, is fully 
determined by the polarization alone, independent of the saturation. The ground-level submatrix p gg is a linear 

combination of the two matrices G and (VW)~ 1 G. For small values of the saturation parameter, the steady-state 
submatrices are 

p ee = SE/a Q , pgg = (V^V^G/ao , (54) 

whereas in the limit of strong saturation S^oowe obtain 

p ee = E/2ai , pgg = d/2ai . (55) 

The matrix for the optical coherence (50) depends on the intensity only through an overall factor (3SI, that is equal 
to £l/a>o in the low-intensity limit, and that approaches zero for strong saturation. 

B. Transitions J g — J — >J e = J + l 

In order to find operators E and G that obey the relations (28) for transitions J g = J — > J e = J + 1, we introduce 
the operators 

L 

V£ b (*)= E (-l) M f^n L _ M (a), (56) 

M=-L 

which arc proportional to the dot product of the spherical harmonic ulm(sl), introduced in (A.l), and the tensor 
Wigner operator 

T? M = E \Ja,Pa)(-iy h -^C^ a Jb _„ b (J b ,» b \ . (57) 

Mo, Mb 

The indices a and b indicate the levels (e or g), and a is a vector in the complex three-dimensional space. We shall 
use the two multiplication rules 



11 



^(efete) = Vlf +1 (e)9r(e) , ^ e (a)^ 2 e | +1 (b) = V^ +1 (b)^(a) . (58) 

The first equation (58) holds for any vector e, and for arbitrary value of the angular momenta J e and J g . The second 
equation (58) is valid for arbitrary complex vectors a and b, provided that J e = J g + 1. These relations follow from 
the multiplication properties of the operators (56) as given in Appendix C, in particular eqs. (C.3) and (C.5). In the 
notation of eq. (56), the coupling operator is given by V = \J (e • e)(2J e + l)/3 V^ es (e). 
In order to identify the operators E and G we introduce the short-hand notations 

V 2 e ] + i(e) = W; V£ +1 (e) = W, (59) 

In addition to V, we introduce the coupling operator 

V = v/(e-e)(2J e + l)/3 Vf e (e) (60) 

Notice that the operators W and V are raising operators, which map substates of the ground level onto excited states. 

The operators W and V are lowering operators. 

Using the relations (58) one can find operators E and G with the properties (28). They are specified by the 
definitions 

E = WW* , G = WW . (61) 
With the notation (59), and the substitution a = e, b = e* in (58), we obtain the identities 

VW = WV , VW = W ] V , (62) 
which indeed prove the first equation (28): 

EV = WW^V = WVW = VWW = VG . (63) 

The second equation (28) is easily verified when the summation is performed over three Cartesian polarization vectors 
ei. By using the second identity (58), with substituted for e, and e for b, one finds 

D\ED q = J2 D\WW^D q = W ]T D q D\W =G (64) 

g=0, ±1 9=0, ±1 9=0, ±1 

where eq. (13) is used in the last step. As indicated in Sec. IIIC, this also proves the commutation rules (21) for the 
steady state in the present class of transitions. 

Now that we have identified the operators E and G with the desired properties, the steady-state density matrix is 
directly obtained in the form (49) and (50). We use the first equation (28), combined with the commutation relation 
(29) for G, and we introduce the operator 

X = {V^Vy^W , (65) 
acting on the states of the ground level. This leads to the expressions 

p ee = psww^ 

p gg = 13 (xX^ + SWW ) (66) 



Alternatively, the operator X is defined by the relation VX = W. The density matrix is indeed diagonal on the 
natural basis, and the commutation rules (21) hold. 

It is illuminating to express the various operators occurring in the equations (66) in terms of basisvectors. The 
commutation rules (29) can be expressed as 
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[WW\ VV^] = , [Ww\v^V}=0. (67) 

The operators V and W depend on the polarization vector e. From the definition (56) one obtains the relation 

fc» 6 (a)] + = {-iy»- J >Vt{**) , (68) 

which shows that, apart from a minus sign, the Hermitian conjugates of W and V are equal to the operators W and 
V with the polarization vector replaced by its complex conjugate e*. Hence, when the polarization vector is taken as 

e* rather than e, the expressions (66) and (65) hold with the replacements V <-> V and W . The commutation 

rules analogous to (67) are then 

[WW,VV]=0, [W^W, VV] = . (69) 

The first commutation rule (67) shows that the natural basis states \(e)j), defined as eigenstates of VV\ are also 
eigenstates of WW\ and the (positive) eigenvalues are called i/?, with i/j positive. In complete analogy to the relation 
(22) between the states \(g)i) and |(e)j) as coupled by V, we notice that the states Wt|(e)i) are eigenstates of the 
ground-level operator W^W, with the same eigenvalue vf. Hence, we denote the normalized eigenstates as \(g)i), so 
that W can be expanded as 

W = J2vi\(e)i){(g)i\. (70) 

i 

From the second commutation rule in (69) it follows that the ground-level states | (g)i) are eigenstates of the operator 

VV , which means that they form the natural basis for the ground level for the polarization e*. From the second 

identity (67) one finds by the same argument that the states W \{g)i) are eigenstates of W W . From the first identity 
(69) it follows that these states form the natural excited basis for the polarization e*, which we indicate as \{e)i). 

Hence, the operators V and W can be explicitly expressed as 

v = Y,Mm&)i\, W = J2»i\(9)i)(Ui\- (71) 

i i 

For symmetry reasons, the values Xi and Vi must be the same as in eqs. (23) and (70). The operator X defined in 
(65) can be expanded as 

X = J2(v i /\ i )\(9)i)((9)i\. (72) 

i 

The steady-state populations of the natural basis states follow from the expressions (66), with the result 

*P=rtS,*P=rt{)f' + S). (73) 

The relations (26) determine the ratio between the population of an excited state of the natural basis set and the 
corresponding ground state. The populations of different excited states of the natural basis set are proportional to 
the eigenvalues i>f. 

In the special case of linear polarization, the polarization vectors e and e* are identical. When we take the 
quantization axis in the polarization direction, the natural basis of states coincide with the Zeeman states \J g ,n) 
which is only coupled to the excited state | J e , fi). The operator W is proportional to the spherical temsor X^f^Q, and 
the corresponding eigenvalues are proportional to the Clebsch-Gordan coefficients cx C}^"|| i ° J _ At . The resulting 
excited-state populations cx |C , j'^^ J _ M | 2 for this special case of linear polarization driving a transition J — > J + 1 
have been indicated before in ref. [15]. 

Calculations of matrix elements of the operators W, W and X are given in the appendix (B) for the natural 
coordinate frame. However, more physical insight can be obtained by expanding the operators X in the form of an 
invariant superposition of the operators V[ 9 
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y (e.e)(2J + 3)^ o 
In order to find the coefficients Cl we use the property (C.2): 

V 1 e °(e)VZ°(e)= J2 E(L,K)V°°(e) 

K=L-1 



J2C L Vi°(e). (74) 



E(L,K) = (-ir^VWLTT){ K j j j L +1 }c 



w°lo ■ (75) 



The recurrent equations for the coefficients Cl follows from the defining relation VX = W 

E(L — 1, L)Cl-i + E{L + 1, L)Cl+i = Sl,2j+i , L = 0, 1, . . . , 2 J + 1 . (76) 

Depending on whether J is an integer or a halfintcger, the odd or even coefficients are equal to zero. In both cases 
the nonzero coefficients Cl are written as 



_ (2L+1K2J + 3) (2J-£)!(2J + L+1)! 

6i "V 3(2j + i) (ITTT)! • (77) 

A remarkable peculiarity of the solution (77), which becomes manifest at large angular momentum J, is a rapid 
decrease of the coefficients Cl with a decrease of index L. This fact allows approximate calculations by restricting 
the expansion (74) to only a few terms. For example, the ratio 



C2J , . , - s / 2 J 



(4J + 1W — - w 2\/2 J 



C27-2 v V 4 J — 3 

for J = 4 (that corresponds to the cycling transition of the D 2 -line of 133 Cs) amounts to about 13.3. As a result, if 
we use the approximation 



then in calculating quantities as the population or the orientation of the levels and the average dipolc moment, the 
error remains below one percent. 

The normalization constants on can be found in explicit form for arbitrary J. From the sum rule for spherical 
harmonics (A. 4) it follows that 

<x x = Tr {wwt} = P 2J+1 (^y) • (78) 

Using the expansion (74), we find 

"° = T ' {«' } " (e.e)(L + 3 ) £*^(<£S))- (79) 

L — 

The coefficients ao and oc\ are simultaneously even (for J a halfinteger) or odd (for J an integer). As one expects, 
the populations of the levels determined by the ratio of «o and cei do not depend on the sign of (e • e). 



C. Transitions J g = J — > J e = J with J a halfinteger 

For a transition between levels with equal half-integer values J of the angular momenta, the coupling operator V 
on the basis of the Zeeman substates is represented by a square matrix. This implies that the concept of the inverse 
operator (V)" 1 as corresponding to the inverse matrix can be used, in the sense that (y) _1 y = n g , ^(V)^ 1 = U e . It 
is nearly trivial to find operators E and G that obey the conditions (28) for this class of transitions. It is immediately 
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obvious that these conditions are satisfied with the choice E — II e and G — H g , in view of the identity (14). Then 
expressions (49) and (50) for the density matrix take the form 

p ee = 0SIl e 

p a9 = P kv^v)- 1 + sn g ] = p livy 1 ^ 1 + sn g ] (so) 



eg 



\' 9e ) {s + ij^y ' 



The commutation rules (29) are trivially obeyed, so that the density matrix is diagonal on the natural bases. These 
expressions have been obtained in Ref. [23] . Remarkably, eq. (80) shows that for this class of transitions, the excited- 
state density matrix p ee is always isotropic at arbitrary field parameters. The ground-state density matrix p gg consists 
of two parts, the term proportional to (V'V) -1 is anisotropic, and describes the distribution among Zeeman substates 
in the low-saturation limit S<1. The isotropic part is dominant in the limit of strong saturation (55). 

In fact, the expressions (80) can also be represented in the form (66), in terms of operators W, W and X. In the 
same spirit as eq. (59), for the present class of transitions we introduce the operators 

V a eB (e) = W: Vf(e) = W, (81) 

where now the rank of the operators takes the minimal value 0. With this definition, we can maintain the expressions 

(61) for the operators E and G, which are indeed isotropic. With the definition (60) of V, the identities (62)-(64) 
remain valid in the present case. This follows from equations (C.3) and (C.4) of Appendix C. With the definition 
(65) of the operator X, the proof of eqs. (66) as well as of the expansions (70)-(72) can be carried over directly. The 
main simplification is that in the present case of a transition J g = J — > J e = J, the coefficients Vi are all the same, 

while the natural basis states \(e)i) and \ (g)i) have an identical vector form. Therefore, on the basis of Zeeman states 

of the excited and the ground level, the raising operators W and W have the form of a unit matrix. 

The matrix elements of the inverse operators (V) -1 , (V^) _1 and (VW) -1 can be found in a closed analytical form 
in the natural coordinate frame for all J (see Appendix B). However, for further applications, it will be convenient to 
express (V) -1 in an invariant form in terms of the spherical harmonics iiiM(e) defined in (A.l). Since the coupling 
operator V is proportional to the operator (56) with the rank L = 1, we can write in the present case 



" \/ (e.e)(2 J+1 ) (for" 



We use the equations 

L+l 



V?[e)V!?(e)= £ E(L,K)V?(e), 



K=L-1 



iKO 
'10 L0 ' 



(82) 



E(L,K) = (-ir+ L +wmL+i){ K j jj}c 

which is a special case of (C.2). The operator (V^ e£, (e)) _1 is expanded in the operators V[ e (e) according to 

(^r^E^ne). (83) 

L=0 

In order to find the coefficients Cl we substitute this expansion into the identity 

FHeX^He))- 1 = n e . 

While noting that E(L, L) = (since Cf ° i0 = 0) one arrives at the two-term recurrent relation 

E{L-\,L)C L - 1 +E(J J + l,L)C L+1 = 5 Lfl y/2J+l, L = 0, 1, . . . , 2 J . (84) 
Expressions for the coefficients with odd indices follow from the relations (84), with the result 
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r ( n^ ^- 1 )" (2^ + 1)2J(2J + 1)(2J + 2) (2J + L)!!(2 J — L — 1)!! 

~ L!^V 3 (2J-L)!!(2J + L + 1)!! (85) 

This is a sequence of terms with alternating signs, which slowly grows in absolute value with the index L. It follows 
from (84) that the coefficients Cl with even index are equal to zero. 
The normalization coefficient (52) is obviously 

ai = Tr{n e } = 2 J + 1 . (86) 

The other normalization coefficient ao = Tr{(V^) _1 } follows from the expansion (83), the orthogonality of the 
Wigner operators and the sum rule for spherical harmonics (A. 4), with the result 

2 7 

«-(^7TT)j:/K<^))- (87) 

This expression contains only even powers of (e • e). 

VI. EXCITED-STATE POPULATION AND AC STARK SHIFT 

In this section, we study the polarization dependence of the total excited-state population. This determines, for 
example, the total fluorescence and thereby also the total absorption in an atomic vapor. For unpolarizcd atoms 
and in the low-saturation limit, the total excited-state population is given by ni unpol) = (5/3) (2 J e + 1)/(2J S + 1) , 
independent of the polarization. This is the case of linear absorption. We now consider as an application of the results 
of the previous section the total excited-state population ir e = Tr{p ee } in the steady state. When we calculate the 
trace of (49), while taking into account normalization (12), we find 

- - Sai/ao (88) 



1 + 2 Sai/ao ' 

where the real functions ot\(e) and ao(e) are expressed in (78) and (79) for the transition J — > J + 1, and by (86) 
and (87) for the transition J — > J (recall that e • e = cos(2e)). As is seen from (88), the analytical expression for the 
total excited-state population is similar as for a two- level atom with non-degenerate levels (see e.g. [1]) but with an 
effective saturation parameter S — Sai/ao, which depends on the class of transitions and on the light polarization. 
If the saturation intensity I sat is defined as the intensity at which n e = 1/4 (at zero detuning) , then we obtain 

_ ao(e) 2tt 2 ftc 7 
lsat{£) - — 7T to , J-0 — ro- , (89) 

where Io is the usual measure for the saturation intensity. This shows that now the excited-state population, and 
thereby the absorption, depends on the polarization. 

This is demonstrated in Fig. 3a. for the class of transitions J — > J with J a half-integer. This shows that the total 
absorption cross-section is reduced compared with the case of linear absorption, so that the medium becomes more 
transparent due to optical pumping. The reduction is the lowest for linear polarization, and it is complete for circular 
polarization. Moreover, the reduction increases with the value of J. Recall that for CPT transitions, J — > J — 1 and 
J — > J with J an integer, the stationary excited-state population is obviously zero, and the transparency is complete. 
The decrease in absorption indicates that the atoms are pumped to states that are more weakly coupled than average. 
Roughly speaking, this implies that pairs of coupled states \ (e)i) and \(g)i) for which the total steady-state population 
is large, tend to have a relatively small coupling constant \\. Since the form (19) of the light shift operators shows 
that the states \{e)i) and \(g)i) are shifted by —SSXf and SSXf, we may also conclude that this class of transitions 
J — ► J with J a half- integer tends to pump the atoms to states with lower AC Stark shift. 

In contrast, for J — » J + 1 transitions the absorption is enhanced by optical pumping. This is shown in Fig. 3b. 
Again, the effect of optical pumping on the total absorption increases with J, and it increases also when the polarization 
is varied from linear to circular. These results are related to the recently discussed effect of clcctromagnetically induced 
absorption (EIA) under two- frequency excitation in the Hanlc configuration on J — > J + 1 transitions [44]. Indeed, 
when the frequencies coincide or the magnetic field is zero, we have a situation close to the stationary interaction 
of atoms with clliptically polarized light considered here. When the frequency difference or the Zeeman splitting is 
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sufficiently large, significant ground-state depolarization appears, and the absorption should be close to the linear 
absorption of unpolarized atoms. In this sense, the enhancement of absorption by optical pumping as illustrated in 
Fig. 3.b resembles EIA. By the same argument as used above, we conclude that J — > J + 1 transitions tend to pump 
the atoms to state with larger AC Stark shifts. This tendency was found before in special cases in the context of 
sub-Dopplcr laser cooling by polarization gradients [11]. 



VII. BROAD-BAND RADIATION 



So far, we discussed the steady-state solutions of the generalized optical Bloch equations (8)-(ll), which describe 
an atomic transition driven by monochromatic polarized light. In this chapter we point out that the results can be 
generalized to the case of light with a finite bandwidth. Broad-band radiation is described by modeling the electric 
field as a stationary stochastic process. The dynamics of an atom in such a field with central frequency w is described 
by the same eqs (8)-(ll), where now the Rabi frequency Q(t) is a complex-valued function of time, proportional to 
the positive-frequency part of the fluctuating electric field. We are interested in the steady-state stochastic average 
of the submatrices p ee and p gg . The time-dependent solution of eq. (8) is 



Peg® 



r 

Jo 



dr exp[(td - 7/2)r]ifi(t - r) V(e)p gg (t - r) - p ee (t - r)V(e) 



(90) 



When we substitute this expression and its analogue for p ge (t) = p eg {t) into eqs. (10) and (11), we arrive at a pair of 
stochastic integro-differential equations for the submatrices p ee (t) and p gg (t). The r.h.s. of these equations contains 
the stochastic parts 

£l*(t)£l(t — T)p gg (t — t), £l*(t)£l(t— r)p ee (t~ r) and their Hermitian conjugates. Stochastic averaging of these equations 
leads to a closed set of equations for the steady-state stochastic averages (p ee ) and (p gg ), provided that the stochastic 
average of these terms may be factorized as 



<fi*(t)fi(t - r)p gg (t - T )) = (JT(t)fi(t - r)) (p gg ) 



(91) 



and similarly for the other terms. When we substitute this factorized form of the type (91) in the equations for p e 
and p gg , we arrive in the steady state at closed equations for the stochastic averages (p ee ) and (p gg ) 



7 (Pee) = -R{VV\ (p ee )} + 2RV (p gg ) t>t + lL [VV\ (p ee )}, 



(92) 



- 7 ]T Dl(% e )D q = -R{V^V,(p g g)} + 2RVU% e )V-tL[V^V,(pgg)}. (93) 

9=0, ±1 

The two real parameters R and L are defined by the equation 

/•OO 

R + iL = / drexp[(W-7/2)r](n*(i)fi(i-r)>. (94) 
Jo 

The quantity R is a measure of the stimulated transition rates, and L determines the strength of the light shift. When 
we substitute this factorized form of the type (91) in the equations for p ee and p gg . 

The factorization (91) is exact in the special case that the finite bandwidth is due to phase fluctuations only. In that 
case we can write Q(t) = Oexp(— iip(t)), with a stochastic phase ip. The factorization is then justified since the phase 
change ip(t) — i/j(t — T) in the time interval [t — T, t] can safely be assumed not to depend on the phase ip at times before 
t — t, which determine p gg (t — r). The phase fluctuations are then described by the independent-increment model 
[45], which has the phase-diffusion model as a special limiting case. The stochastic average of the field correlation 
function then decays exponentially, according to the equality (n*(t)n(t — t)) = |f2| exp(— fir), with p the bandwidth 
(halfwidth at half maximum) of the Lorentzian profile. In this case the quantities R and L are determined by the 
equation 

R + iL= . (95) 

p + 7/2 — iS 

It is a simple check to notice that eqs. (92) and (93) reduce to the corresponding equations (16) and (17) for 
monochromatic light when we substitute p — 0. In that case we simply find R = "fS/2 and L = SS. 
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In the case that the driving light also has intensity fluctuations, the situation is more complex, and the factorization 
(91) is not exact. When the fluctuations are sufficiently weak and sufficiently rapid, we can still assume the factoriza- 
tion as a reasonable approximation. Therefore, the equations (92) and (93) can be assumed to be valid for broadband 
radiation in many situations of practical interest. These equations, which strongly resemble the corresponding equa- 
tions (16) and (17) for monochromatic light, determine the steady-state stochastic average of the density matrices p ee 
and p gg . In the preceding sections we have demonstrated that for monochromatic light and for all allowed values of 
J e and J g the steady-state density matrix obeys the commutations rules (21), so that the density matrix is diagonal 
in the eigenstates of the light-shift operators. This means that the solutions of (16) and (17) do not depend at all on 
the strength of the last terms in these equations. We conclude that the expressions for the stochastically averaged 
steady-state solutions (p ee ) and (p gg ) in the absence of a dark state coincide with the solutions obtained in Sec. V 
for p ee and p gg , with the simple replacement S — > 2R/-f. The structure of the dark states follows from the defining 
equation (31), so that also the results of Sec. IV are not modified by the fluctuations of the driving light. The 
steady-state polarization properties of the atom are basically unaffected by the finite bandwidth. The steady-state 
optical coherences p eg are best described by the expression for the stochastic average 

(fl*p eg ) = (iR - L) [9(e) (p gg ) - (pee) V(e) 

which follows immediately from eq. (90). An expression containing p ge follows after Hermitian conjugation. Obviously, 
these conclusions are valid exclusively when the light polarization displays no fluctuations. 



VIII. DISCUSSION AND CONCLUSIONS 



We have given a complete analytical and invariant description of the steady state density matrix of a closed atomic 
dipolc transition J g — > J e driven by a resonant polarized radiation field. This is a long-standing problem in atomic 
and optical physics. Solutions have been known for some time in special cases of polarization and values of the 
angular momenta J e and J g of the excited and the ground state. The most complex class of transitions occurs for 
Je = Jg + 1- I n this case, the excited-state density matrix can be highly non-anisotropic. It is remarkable, however, 
that the anisotropy depends exclusively on the ellipticity of the polarization, and it is unaffected by the frequency 
detuning of the radiation from resonance, the light intensity or the spontaneous decay rate. In the case that J e = J g 
is half-integer, the excited-state density matrix is fully isotropic in the steady state. In the remaining classes of 
transitions (J e = J g is integer, and J e = J g — 1), the system has one or two dark states, and the degree of excitation 
vanishes in the steady state. For these cases, we give analytical invariant expressions for these dark states for arbitrary 
elliptical polarization. These results arc interesting, not only from a fundamental point of view as an exact solution 
of a quantum mechanical problem, but also since they can be used in numerous applications. 

As a first example, we mention the problem of non-linear propagation of elliptically polarized light in a resonant 
gas medium. The steady-state solution allows one to find the non-linear susceptibility tensor in analytical form. The 
Doppler broadening is taken into account by the substitution 5 — > 5 — k- v in the expressions for p eg , and then average 
over velocity. 

A second case of interest is high-resolution polarization spectroscopy. The Doppler-free resonances in the scheme of 
a strong pump and a weak probe field can be directly evaluated by calculating the linear response to the probe, in a 
steady state that is determined by the pump. Non-linear interference effects between pump and probe are negligible 
in several cases, e.g. when they are counter-propagating. 

A third situation of practical importance occurs when cold atoms are slowly moving through non-uniform radiation 
fields, with a position-dependent amplitude Eo(r) and polarization vector e(r). The steady-state solution discussed 
in this paper can be viewed as the zeroth-order approximation with respect to the atomic velocity. This solution is 
needed for an explicit calculation of radiative forces [19,22,41,46] and geometrical potentials [42,47], which also affect 
the dynamics of atoms in optical lattices. 

Generally speaking, in many problems there are factors not taken into account in our solution, such as finite 
interaction time, translational motion of atoms, magnetic field etc. Very often these factors can be considered as 
a small perturbation. In all these cases the steady-state solution presented in this paper constitutes a zeroth-order 
approximation, and thereby the first necessary step in the corresponding perturbation treatment. 
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APPENDIX A: SPHERICAL HARMONICS OF A COMPLEX DIRECTION 

Throughout the paper we use spherical harmonics that differ from the standard definition [35] by a multiplicative 
factor 



4ir 

nLM = xl 2LTi YLM ■ 

For arbitrary complex vector a = a' + ia" a spherical harmonic of the rank L is defined in terms of the tensor 
constructions (33): 



n LM (a)^y (2J ^ 1)!! {a} L , (A.l) 

where a = ^/ (a • a). These generalized spherical harmonics depend only on a direction in the complex three- 
dimensional space [37], i.e. they do not change under the transformation a — > va with v an arbitrary complex 
number. For real vectors (a" = 0) the definition (A.l) leads to the standard spherical harmonics [35]. Starting from 
equation (A.l), one can derive the well-known formula [35]: 



n LM (a) = e ^y||_^£||pf ( cos0) , (A.2) 

where Ppf(x) are the associated Legendre functions, and the complex parameters 9 and <p ar e expressed in terms of 
the spherical components of vector a by 

cos6> = ao/a , e 2 ^ = — a + i/a_i . 

Formula (A.2) can be regarded as a suitable analytic continuation of the standard definition of the spherical harmonics 
nijf(^) [35] to complex values of the angles </> and [38]. The definition (A.l) is important since the functions 
n-LM obey the same group-theoretical relations as ordinary spherical harmonics [38]. In particular, we indicate the 
Clebsch-Gordan expansion of the product of two spherical harmonics of the same argument 

n hmi (a)n hm2 (a) = J2 C h0 ho^hmr bm^iMfa) (A.3) 

LM 

and the sum rule for the dot product of spherical harmonics of different arguments [38,37] 

(n L ( a )-n L (b))=P L ( { ^P-) , (A.4) 



ab 



where Pl{x) are the Legendre polynomials. 



APPENDIX B: CALCULATING MATRIX ELEMENTS 

The matrix elements of operators V -1 , i W> and X, which are used to write the steady-state density 

matrix p, can be determined in the natural coordinate frame. To be specific, we fix the sign in cq. (3): 

e = ^/cos(2e)e — \/2sin(e)e+i. (B.l) 
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1. Transitions J g — J — > J e — J with J halflnteger 

In the natural coordinate frame the matrix V is real and has a lower triangular form with two nonzero diagonals: 

/ • 0\ 



V 



VfJ,,fJ,— l */J>tJ. 



where in accordance with the definitions (7) and (B.l) 



(B.2) 



v, 



l (J + fl)(J-y + l) 
J(J+1) 



sin(e). 



(B.3) 
(B.4) 



Its inverse matrix also is of the lower-triangular form and real. The matrix elements of V 1 are calculated by a direct 
method: 



V~ 



(-iy 



n 



(a-1) 



<J(J+1) / sin(e) 



cos 



j_ A y/(J + a)(J-a + l) 

(2 £ ) vv^y; M' Q= y +1 



(B.5) 



The repeated products in (B.5) should be read while using the conventions 



n fa = !; n f° = ° if > 

Since the matrix V is real, is obtained from V^ 1 by transposition, i.e. 

can easily write the matrix elements of : 



(vt) 



. Thus, one 



y a {a-\) \ ( T~T K*'(a'-1) 



n 



E 



sin(e) \ 1 



J(J+1) 

cos(2e) \ ^cos(2 £ ) ^ ^ 2 



n 

Va=i/+1 



V(J + a)(J-a+l) 



a 



^(J + a'XJ-a' -;- 'IT 

i a'— z^+1 



a' 



(B.6) 



2. Transitions J g = J — > J e = J + 1 

In the natural coordinate frame the components of the spherical harmonics (A.l) of the polarization vector e are 
written as 
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if M > 0, and ul-m (e) = for M < 0. Substituting (B.7) into the definition (56), we arrive at (J a = J + 1, Jb = 
J, L = 2J+ 1): 



j_ ro (2J+l + /i-m)l / 

V 



(2J + 2)!(2J)! 



(4J + 1)!(J + 1 + M )!(J + 1 - /i)!( J + m)!(J - m)! 



\ fl — 771 

' ] , (Bi 



\ ^/cos(2e) / 

where ^ = — J — 1, — J, ... J + 1, m = — J, — J + 1, ... J and ^ — to > 0. The matrix IT can be obtained from (B.8) 
using the time-reversal operation W mfl = (—l) J s~ m ~ Je ~ fl W_ / j j _ rn . 

In order to find matrix elements of X we decompose it X = UW using the Moore-Penrose pseudoinverse [39] matrix 
U with respect to V, i.e. VUV = V. The nonzero elements of V are given by 



V ^ = V (jTWjTT) Vcos{2£) 



^- 1 = V (J+1)(2J+1) bm£ - (R9) 



As the pseudoinverse matrix to V we take the matrix with elements 



TJ I (J+1K2J+1) / sine \ m ~" ™ / J + ^ , . 

mAI ~V( J + 1 +^)( J + 1 -M)co S (2 e ) ^ yc^iyj Ji^j+i-i/' (ti - W) 

for /x = — J, — J+ 1, ... J, supplemented by the zero columns t/ m _ j_i = and U m j + \ = 0. The matrix multiplication 
of (B.10) by (B.8) yields the final result for the elements of X: 



/(J + 1)(2J+1)!(2J + 2)!(J + to)!(J-to)! / sins \ 

X m m> = \ — TTTT^ ^TT^ — I , I 



m 



(4J+1)!(J + to')!(J-to')!cos(2 £ ) ^ A /cos(2e) < 

r Vl I-'--" (2J+l + /i-TO')! 

^-, ( lj ( M -m')!(J+l + )u)!(J+l-//)! • (RU) 



APPENDIX C: ALGEBRA OF THE OPERATORS V^ B (A) 

Using the standard Racah algebra [40], one can write a general expression for products of the operators V (56) with 
different ranks: 

% fc (a)^(b)=^(-l)^+^n Ll , L2 |^ ^ j 2 a }({n Ll ( a )®n L2 (b)} K -T£) , (C.l) 



where Tl x , y ,... — \J (2a; + l)(2y + 1) . . . is the standard notation of the handbook [35]. In the special case that b = a, 
after using eq. (A. 3) we obtain from (C.l) an analogue of the Clebsch-Gordan expansion: 

%"(a)%(a) = EC" 1 ) J " +Jc+Ll+L2 n i ,L 2 Cf 1 ° L20 { \ L l ^ } Vn*) • (C.2) 
Equations (C.l) and (C.2) lead to the following relationships: 

1) For arbitrary ranks L\ and L 2 , and for arbitrary angular momenta J a and Jb we find 
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V£(a)V£(a) = V£(a)V£(a) , (C.3) 

since both sides have the same expansion in the tensor operators V^ a (a). Here we use the symmetry of (C.2) with 
respect to the permutation L\ «-> L2 at J a = J c . 

2) Depending on the class of transition, for arbitrary vectors a and b we find: 

a) for transitions J g = J — > J e = J 

Vr(h)V Q ea ( a ) = ^ e (a)^(b) . (C.4) 

b) for transitions J g = J — > J e = J + 1 

^ e (a)^f +1 (b) = i>/; +1 (b)^(a) . (C.5) 

c) for transitions J g = J — > J e = J — 1 

?i efl (a)V&-i(b) - ^|- 1 (b)^r(a) . (C.6) 

The property (C.4) is obvious, if we recall that in this case the operator Vq C is proportional to the unit matrix and 
Vf e (h) = Vi 9 (h). To prove the validity of (C.5) it is sufficient to expand both sides of (C.5) in the operators Tf? q 
and allow for the fact that all ranks except K = 2 J are forbidden by the selection rules contained in the 6j symbols 
in (C.l). The equation (C.5) then reduces to the identity {ni(a) ® ri2j+i(b)}2j = {ri2j+i(b) ® ni(a)}2j, which holds 
since the number (2 J + 1) + 1 — 2 J = 2 is even. The equation (C.6) can be proved in a similar way. 
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FIG. 2. Natural coordinate frame for the representation of elliptical polarization, a) polarization vector is superposition of 
circular polarization denned by cylinder (in x'y'-plane), and linear polarization along axis of cylinder (V-axis). b) transition 
scheme with the z'-axis as quantization axis. 
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FIG. 3. Total absorption at low saturation versus the ellipticity e. (a) transitions J — ► J with J a half-integer; J runs 
from 1/2 to 9/2 (from top to bottom), (b) transitions J — > J + 1; J = 0, 1/2, 1, 2, 3, 4 (from bottom to top). All curves are 
normalized to the linear absorption. 
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